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Abstract 

Microscopic optical model potential results for reaction cross sections of proton elastic scattering 
are presented. The applications cover the 10-1000 MeV energy range and consider both stable 
and unstable nuclei. The study is based on in-medium g-matrix full-folding optical model 
approach with the appropriate relativistic kinematic corrections needed for the higher energy 
applications. The effective interactions are based on realistic NN potentials supplemented with 
a separable non-Hermitian term to allow optimum agreement with current NN phase-shift 
analyzes, particularly the inelasticities above pion production threshold. The target ground-state 
densities are obtained from Hartree-Fock-Bogoliubov calculations based on the finite range, 
density dependent Gogny force. The evaluated reaction cross sections for proton scattering are 
compared with measurements and their systematics is analyzed. A simple function of the total 
cross sections in terms of the atomic mass number is observed at high energies. At low energies, 
however, discrepancies with the available data are observed, being more pronounced in the lighter 
systems. 

PACS numbers: 25.40. Cm, 25.60.Bx, 24.10.Ht 

Keywords: Proton scattering, optical model, total reaction cross-section 



*Electronic address: arellano@dfi.uchile.cl; URL: 



http : / / www . omp- onl ine . c 1 



I. INTRODUCTION 



Nucleon-nucleus integrated cross sections constitute a key observable in fundamental nu- 
clear research as well as in applications of nucleon-induced reactions. These quantities are 
of particular importance in nuclear technology such as nuclear transmutation, nuclear waste 
treatment, safety assessment and medical therapy, among many examples. From a funda- 
mental point of view, their description within global microscopic approaches constitutes a 
stringent test to the understanding of the underlying physics involved in the interaction of 
nucleons colliding with a nucleus. Depending on the energy range of the projectile, total 
cross sections may constitute an important input in the study of diverse phenomena. For in- 
stance, the study of r-processes in astrophysics require the knowledge of total cross sections 
at energies near and below a few MeV, whereas spallation applications may well require data 
up to GeV energies. 

In this article we present a study of the systematics exhibited by proton-nucleus (pA) 
reaction cross sections based on the microscopic full-folding (FF) optical model potential 
approach The study spans in energy from a few MeV up to 1 GeV, considering 

various even-even isotopes from carbon up to lead. The calculated reaction cross sections 
are compared with existing data, whereas results for proton scattering from unstable nuclei 
are examined as functions of the number of constituents of the target. 

Research on proton reaction cross section of pA elastic scattering has received significant 
attention during recent years. These studies include experimental 3J, phenomenological 
0, Q and microscopic 6, 7, ^, q| analyses. Global analyses in the framework of Glauber 
;heory have also been reported 10|. Recent experimental efforts by Auce and collaborators 
3I have been very valuable in reporting new measurements at intermediate proton ener- 
gies. The phenomenological studies of Refs. {4], Q represent attempts to provide simple 

geometry and mass distribution 



3, |8| are based on the referred 



parametrizations of the observed cross sections in terms of 
of the targets. The microscopic studies reported in Refs. 6 
'coordinate-space g-folding model' [ill ], where a local medium- dependent NN effective inter- 
action is folded with the target ground-state mixed density. As a result, a nonlocal optical 
potential is obtained, where the non locality stems from the inclusion of the exchange term. 
The results reported here differ from the above g-folding model in two aspects. The first of 
them is that the antisymmetrized NN effective interaction is handled in momentum space 
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with no special considerations on an eventual local structure in coordinate space [lj. Thus, 
we retain the intrinsic off-shell behavior of the interaction. The second difference lies in 
the representation on the mixed density. Whereas the ^-folding approach treats explicitly 
the mixed density extracted from shell-models, the results we report here rely on its Slater 
representation using only the diagonal elements. An assessment of this approximation has 
been made within the free t matrix full-folding approach, where it is shown that it affects 
only slightly the differential scattering observables at momentum transfers (q) of 1 fm _1 , 
becoming more visible at q > 2.5 fm _1 (c.f. Fig. 3 of Ref. 12j). 

This article is organized as follows: In Section II we outline the general framework upon 
which we base our study. In Section III we present results for the calculated reaction cross 
sections for pA elastic scattering at energies between 10 MeV and 1 GeV. Additionally, we 
examine the systematics exhibited by the calculated reaction cross sections and the total 
cross sections for neutron- nucleus (nA) elastic scattering at energies above 400 MeV. Finally, 
in Section IV we present a summary and the main conclusions of this work. 



II. FRAMEWORK 

The microscopic approach we follow is the in-medium FF optical model approach [lj], 
which is a realization of the double convolution of an NN antisymmetrized effective inter- 
action with the target ground-state mixed density. In its formulation, with the use of the 
Slater representation of the mixed density and the assumption of a weak dependence of the 
g matrix on one of the momentum integrals, the FF optical potential for pA collisions at 
beam energy E takes the simplified form 

U pp (k',k;E) = ! dZ e iik '- k) z ^ x 

a=p,n 

p a (Z)g pa [k',k;p A (Z)} , (1) 

where p a is the point density of specie a, and g pa represents an off-shell Fermi-averaged 
amplitude in the appropriate pa pair and evaluated at the local target density pa{Z). More 
explicitly, in an infinite nuclear matter model these density-dependent amplitudes are given 
by 

g NN (k',k;p) = -X- [ Q(k-\P\)g(k r ',k r ; ^s; p) dP , (2) 



3 



where g(k r ', k r ; \fs\ p) corresponds to off-shell g matrix elements for symmetric nuclear mat- 
ter of density p. The relative momenta, k r and k r , depend on the asymptotic momenta 
k', k as well as P, the mean momentum of the struck target nucleon. Their general form 
is k r = Wk-(l- W)p, with p = P/2 + k - k', and W = W(E; k', k, P). As there is no 
local prescription to handle the g matrix, the results we report here retain the intrinsic non 
localities of the NN interaction. The s-invariant specifies the energy at which the interaction 
is evaluated, which depends on the beam energy E and kinematics of the interacting pair. 
Details about the implementation of the Fermi-averaged elements and relativistic kinemat- 



ics corrections have been discussed in Ref. 



• 3- 



With the above considerations we obtain a 



nonlocal optical potential which is handled exactly within numerical accuracy. 

An important element for the realization of the optical potential is the two-body effective 
interaction, which in our model takes the form of the nuclear-matter g matrix. This is based 
on the Brueckner-Bethe-Goldstone model for symmetric nuclear matter and described by 
the integral equation 

g(u) = V + V ? — - g(u) , (3) 

uj - e x - e 2 

with V the free-space two-nucleon potential, e\ and e 2 the quasiparticle energies, Q the 
Pauli blocking operator and u the starting energy just above the real axis. In the case of 
the free t matrix applications, the Pauli blocking operator is set to unity and the nuclear 
self-consistent fields vanish. Thus, g — > t, the NN scattering matrix in free space. 

A physically acceptable g matrix for GeV nucleon energy applications requires a bare AW 
potential model consistent with its high energy phenomenology, particularly the presence of 



complex phase-shifts above pion-production threshold 



131 ] . In our approach this is achieved 



by supplementing realistic NN potential models (Vr) with a separable term of the form 

v = v R +\t)r(E)(z\ . (4) 

In this study the form factors | £) are taken as harmonic oscillator states characterized with 
ftw =450 MeV. The strength T(E) is energy and state dependent, becoming complex in 
those states where loss of flux is observed. These coefficients are calculated analytically for 
each state to reproduce the year-2000 np continuous-energy solution (SP00) of the phase- 
shift analysis by R. Arndt et al. {3]. With these considerations we are able to account for 
the absorption stemming from the elementary AW interactions. Since most realistic bare 
potentials do not share a common phase-shift data basis, in particular the SP00 solution, we 



suppress the separable contribution below 300 MeV. The inclusion of the separable strength 
is done gradually out to 400 MeV, energy from which the full strength is taken into account. 

Another important input in our calculations is the target ground-state radial density. 
These are obtained following self-consistent Hartree-Fock-Bogoliubov calculations with the 
Gogny force Q Q. Th,s interaction contains a centra. finite-range tern, a grange 
spin-orbit contribution and a zero-range density- dependent contribution. This approach has 
been instrumental for obtaining the spherical radial densities for 12 C and the isotope families 

12-26 0; 34-64 Ca) 50-86^ 80-100^ 96-136g n and 176-224p b induded in this study . 

III. APPLICATIONS 

A. Scattering from stable nuclei 

Following the considerations outlined above, we have calculated the reaction cross sections 
for proton-nucleus elastic scattering an at beam energies ranging from 5 MeV up to 1 GeV. 
The targets considered are 208 Pb, 90 Zr, 60 Ni, 58 Ni, 40 Ca, 16 and 12 C, whose proton and 
neutron root-mean-square (r.m.s.) radii are summarized in Table [H The results for an are 
shown in semi- log scale in Fig. [TJ where the data are taken from Refs. jsj (open squares), [3] 
(filled circles) and (triangles). The calcium and nickel data have been colored blue and 
red, respectively, to distinguish them when they overlap. The solid curves represent results 
from the FF optical model using the nuclear matter g matrix based on the np Argonne V18 



reference potential [18], whereas the dashed curves represent results using the free t matrix 
corresponding to the same bare interaction. As observed, the level of agreement between 
the calculated cross sections and the data is qualitatively different above 200 MeV from 
that below 100 MeV. At proton energies below 100 MeV -with the exception of 208 Pb- all 
FF results overestimate or, being more pronounced for the 12 C and 16 targets. In these 
two cases differences of ~100 mb at 30 MeV are observed, with no proper account of the 
maxima near 25 MeV exhibited by the data. Instead, the calculated cross sections decrease 



monotonically in the range 10-300 MeV. This feature is also observed in Ref. [9(, where a 
much closer agreement with the data is reported. The deficiencies of microscopic optical 
models at low energies have been addressed recently in Ref. [l9|], where it is suggested that 
the inclusion of nucleon-phonon couplings may be needed to improve the agreement with 
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the low energy data. 

For the heavier targets ( 40 Ca, 58 > 60 Ni, 90 Zr and 208 Pb) the calculated cross sections follow 
the same qualitative trend of the data. The discrepancies observed at energies between 20 
and 100 MeV can be characterized by a uniform 50-mb overestimate of the calculated cross 
sections relative to the data. At energies above 200 MeV the agreement with the measured 
cross sections is considerably improved, with the exception of the 860-MeV measurements 
of 12 C and 208 Pb, where a clear disagreement is exposed. At this specific energy we have 
performed FF optical model calculations for the targets considered and our predictions are 
summarized in Table IHII These predictions are consistent with the data from Ref. near 
and above 1 GeV (open triangles), an indication of the adequacy of our approach at these 
high energies. 

For completeness in this section, in Fig. [2] we compare the measured {20] and calculated 
nA total cross sections (cry). The calculated results were obtained within the FF optical 
model approach, and the applications cover the energy range between 5 MeV and 1 GeV. 
What becomes clear from this figure, in contrast to the description of <tr, is that the clos- 
est agreement with the data occurs with the lighter self-conjugate targets. Instead, for the 
isospin asymmetric 90 Zr and 208 Pb systems, the calculated <7t lacks the pronounced oscillat- 
ing pattern exhibited by the data. Microscopic calculations within the ^-folding approach 
have been reported |9( to provide an excellent account for the total cross-section data in the 
60-200 MeV energy range. Its extension to 600 MeV has been achieved by means of a simple 
parametric form 2jJ . From the prospective of the microscopic FF approach, the account for 
or for isospin asymmetric targets remains a pending issue. 



B. Scattering from unstable isotopes 

In order to explore the behavior of the integrated cross sections as functions of the mass 
number, particularly its isotopic asymmetry, we have evaluated or and <jt for nucleon- 
nucleus elastic scattering. These applications include results at 0.4, 0.7 and 1.0 GeV nucleon 
energies, covering even-even isotopes of oxygen, calcium, nickel, zirconium, tin and lead. In 
Fig. [3] we show a log-log plot of the calculated an (red circles) and ot (blue circles) at 0.4, 
0.7 and 1.0 GeV nucleon energy as functions of A, the number of nucleons of the target. 
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The straight lines represent the least-square regression of the type 

o = a A p , 

with o"o the reduced cross section and p the power-law exponent. Both parameters depend 
on the nucleon energy E. For clarity in the figure, the results for 0.7 GeV and 1.0 GeV have 
been off-set by factors of 10 and 100, respectively. 

As observed from the figure, with the exception of gr for the oxygen isotopes, the overall 
trend of the calculated cross sections follows very well the AP power law. In the case of the 
oxygen isotopes at 1.0 GeV, a slight deviation as a function of the asymmetry A — 2Z is 
observed. In this particular case the maximum deviation of the calculated or with respect 
to the A p law is bound by 3%. Our estimate is that this manifestation of the isotopic 
asymmetry is genuine, although a more reliable evaluation would require a better handling 
of the target mixed density. Such considerations go beyond the scope of the present work. 

In Table Hi] we summarize the results from the power law regression of the calculated 
cross sections. In all cases we include the standard deviation of the obtained reduced cross 
section and corresponding exponent. Judging by the standard deviation, the results at 1 
GeV exhibit the closest agreement with the A p behavior, case in which gr ~ ^0.644^ an( j 
g t ~ A 0,748 . Clearly these results differ from the A 2 / 3 geometric law, suggesting the relevance 
of the hadron dynamics and implicit correlations in the collision phenomena. 

A closer comparison between the calculated cross section and the least-square power 
law fit is shown in Fig. HJ where the solid curves represent the cross sections for nucleon 
scattering from 16 at 1 GeV using the corresponding parameters from Table [Til The filled 
circles represent the calculated gr and gt using the FF approach, where dashed lines are 
drawn to guide the eye. In this figure it becomes evident the departure from the A p behavior 
when the neutron number varies. Indeed, the calculated cross sections are weaker than those 
obtained from the power law for increasing neutron excess. Conversely, the oxygen isotopes 
with fewer neutrons than protons yield cross sections greater than the ones prescribed by the 
parametrization. These features become more pronounced in gr than in gt- Although we 
do not have thorough interpretation of this feature, we notice that the elementary total cross 
sections, g pp and G pn , grow very rapidly between 500 MeV and 1.3 GeV nucleon laboratory 
energy. This increase is more pronounced in a pp than in a pn , where near 1 GeV a pn is 
overtaken by o pp . The smaller a pn relative to a pp weakens the increase of the total cross 
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section as the number of neutrons is increased. 



IV. SUMMARY AND CONCLUSIONS 

We have presented a global study of the reaction cross section for proton scattering from 
various unstable and stable isotopes at energies between 10 MeV and 1 GeV. The study is 
based on the microscopic in-medium FF optical model potential. The effective interaction 
is represented by the nuclear matter g matrix which is obtained from bare AW potentials 
with complete account for the loss of flux of the interaction above pion production threshold. 
This feature is implemented with the inclusion of a separable, energy-dependent component 
added to the reference potential, in this case the Argonne AV18 NN potential model. To this 
purpose, we have used the SP00 phase-shift analysis available from the George Washington 



University, Data Analysis Center [13j. The target ground-state densities were obtained from 
Hartree-Fock-Bogoliubov calculations based on the finite range, density dependent Gogny 
force. 

The calculated (Jr for proton scattering from stable nuclei are in reasonable description 
of the data above 200 MeV, but lack comparable agreement with the data below 100 MeV. 
These discrepancies become more pronounced in the case of the lightest targets studied, i.e. 
12 C and 16 O, but diminish with increasing size. At higher energies the FF optical model based 
on the described NN interactions are consistent with the measured data. A comparison of 
the FF optical model approach with the total cross section data for nA scattering shows 
close agreement in the cases of self-conjugate targets, but clear disagreement for the isospin 
asymmetric nuclei. Thus far we have been unable to identify a microscpic mechanism able 
to account for such discrepancy. 

We have also investigated the systematics of the total cross sections as a function of the 
mass number for various isotope families. The results for both <jr and <jt for proton and 
neutron scattering, respectively, suggest an A p power-law accurate within 3%. At nucleon 
energies between 0.7 and 1.0 GeV the exponents for <tr and <jt are close to 2/3 and 3/4, 
respectively. Some slight deviation from this law is observed when the number of neutrons 
departs from that for the most stable isotope. This feature becomes more pronounced in 
the case of or at 1 GeV for 16 O, but weakens for heavier targets and at lower energies. 

The work reported here constitutes a global assessment of the microscopic momentum- 
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space FF optical model approach on its account for proton reaction and neutron total cross 
section, covering nearly three orders of magnitude in energy. Similar studies have been 
reported within alternative theoretical approaches, such as Glauber theory [h]], the in- 
medium (^-folding approach and global Dirac phenomenology [7| . Although some differences 
occur in the quality of the description of the data, particularly with respect to the q-folding 
approach, it remains difficult to identify the sources of such differences. Indeed, in Ref. [11 1 
the effective interaction is represented in coordinate space as an expansion of Yukawa form 
factors of various ranges and complex energy dependent strength, in the form [c.f. Eqs. (7.1) 
and (7.2) of Ref. 



with 



(i)r 



tfjfM^sf^,— . 

The calculations reported here do not make any consideration regarding the coordinate 
space structure of the g matrix. Actually, they are taken directly from the solution of the 
Brueckner-Bethe-Goldstone equation. Additionally, the representation of the mixed density 
may be critical in the case of the light systems. In this regard, the ^-folding approach is more 
detailed by making explicit use of the full mixed density based on shell models, a treatment 
which is pending in the FF approach. 

From a broader prospective, the realization of the FF optical model relies on a weak 
dependence on the struck nucleon momentum. This assumption is crucial to reduce the 
number of multidimensional integrals by three, making computationally feasible the eval- 
uation of optical potentials in the FF and g-folding approaches. However, the conditions 
under which this assumption becomes most (or least) adequate have not been investigated. 
Other aspects such as charge-symmetry of the bare AW interaction and asymmetric nuclear 
matter effects may also need to be examined. 
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FIG. 1: (Color online) Measured 



13] and calculated reaction cross section based on full- 



folding optical model potential using the g matrix (solid curves) and the free t matrix (dashed 



curves . 



TABLE I: Root-mean-square radii of the point proton (R p ) and neutron (R n ) densities used in this 
study. 



Nucleus 


208 pk 


90 Zr 


60 M 


58 M 


40 Ca 


16Q 


12 C 


R P [fin] 


5.437 


4.210 


3.717 


3.695 


3.405 


2.676 


2.419 


Rn [fm] 


5.573 


4.265 


3.738 


3.683 


3.365 


2.656 


2.402 
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TABLE II: Results for the AP power law regression of the calculated cross sections as a function 
of the nucleon energy. 





OR 


<7y 


E [GeV] 


ao [mb] 


P 


ct [mb] 


P 


0.4 


42.0 ± 0.1 


0.692 ± 0.003 


40.9 ± 1.1 


0.789 ± 0.005 


0.7 


57.0 ± 0.8 


0.652 ± 0.002 


61.4 ± 0.8 


0.746 ± 0.002 


1.0 


59.8 ± 0.1 


0.644 ± 0.002 


63.0 ± 0.1 


0.747 ± 0.001 
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FIG. 3: (Color online) Calculated ctr and gt for nucleon elastic scattering from O, Ca, Ni, Zr, Sn 
and Pb even-even isotopes. The lower-, middle- and upper-pair curves correspond to 0.4, 0.7 and 
1.0 GeV nucleon energy, respectively. 
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TABLE III: Measured [iff] and predicted or for pA scattering at 860 MeV. 



Target 


Measured ctr [ mb ] 


Predicted ctr [ mb ] 


208p] 3 


1680 ± 40 


1852 


90 Zr 




1079 


eo Ni 




823 


58 Ni 




804 


40 Ca 




647 


16Q 




357 


12 C 


209 ± 22 


284 
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